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ABSTRACT 



We investigate the possibility that mass transfer early in the evolution of a massive 
binary can effect a reversal of the end states of the two components, resulting in a neu- 
tron star which forms before a black hole. In this sense, such systems would comprise 
the high-mass analogues of white dwarf-neutron star systems such as PSR B2303+46. 
One consequence of this reversal is that a second episode of mass transfer from the 
black hole progenitor star can recycle the nascent neutron star, extending the life of 
the pulsar. 

An estimate of the formation rate through this channel is first performed via a 
simple analytical approach, and then refined using the results of the SeBa binary 
evolution package. The central role of kicks in determining the survival rate of these 
binaries is clearly demonstrated. The final result is expressed in terms of the number of 
field pulsars one can expect for every single neutron star-black hole (ns,bh) binary. We 
also calculate this figure for black hole-neutron star (bh,ns) systems formed through 
the usual channel. Assuming kicks drawn from the distribution of Cordes & Chernoff 
(1998), we find an expectation value of one (ns,bh) binary per 4 x 10 4 pulsars, and one 
(bh,ns) system for every 1500 pulsars. This helps to explain why neither system has 
been seen to date, though it suggests that detection of a (bh,ns) binary is imminent. 
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1 INTRODUCTION 

The "recycling" of pulsars through mass accretion from a 
binary companion has be en established in t he literature for 
more than two decades dAfpar et al.lll982l) . In the case of 
neutron star-neutron star binaries, the recycled pulsar is in 
fact the first neutron star to form, as the mass accretion is 
driven by Roche-lobe overflow during the evolution of the 
second, initially less massive star. There are several double 
neutron star systems (ns,ns) currently k nown, which have all 
experienced some degree of recycling iArzoumanian et alJ 
ll999l:lBurgav et al.ll20od) . 

Th e second star need not end up as a ne utron star, of 
course . Ivan Kerkwiik fc KulkarnH il999T) and iDavies et alJ 
J2002F) showed that the binary PSR B2303+46, which 
contains a slowly rotating pulsar, was consistent with a 
white dwarf-neutron star (wd,ns) system where the white 
dwarf evolved first. This reversal leaves the neutron star 
with no source of matter to accrete, hence no recycling 
takes place in this instance, and one is left with a white 
dwarf and an unrecycled pulsar. The mechanism behind 



this reversal is detailed in iTutukov fc YungelsorJ Jl993>) 
iPorte gies Zwart fc YungelsorJ <ll999T ) and ^mri^£ Sennelsl 
(2000). Essentially, as the primary evolves and fills its Roche 
lobe, it spills so much mass onto the lighter companion that 
the evolutionary destinies of the two stars are reversed. The 
initially massive primary has insufficient mass to become a 
neutron star, ending up as a heavy white dwarf, whereas 
the companion gains enough mass to evolve into a neutron 
star. We label such a binary as a white dwarf-neutron star 
(wd,ns) system to distinguish its heritage from the more 
common (ns,wd) binaries, in which the neutron star forms 
first. For the (wd,ns) evolution channel, the stars have to be 
of similar mass, and these masses have to straddle or lie just 
below the boundary between white dwarf and neutron star 
formation. 

An analogous situation exists in higher mass systems, 
where both components have a ZAMS mass below the 
threshold for black hole formation. Unlike the (wd,ns) sys- 
tems above, the pulsar in this scenario will experience recy- 
cling if the system undergoes Roche-lobe overflow when the 
companion leaves the main sequence. As yet, no neutron 
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stars have been discovered with black hole companions. An 
estimate of the rate at which such objects are formed, cou- 
pled with educated guesses about the longevity of such a 
system, allows us to infer the number of such systems that 
may exist in our Galaxy. A more useful comparison is the 
ratio of the number of (ns,bh) systems to pulsars, as an 
expectation value for the (ns,bh) detection rate per pulsar 
found. If this number is greater than the number of pulsars 
currently known, the absence of (ns,bh) binaries is easily un- 
derstood. Conversely, if the rate implies that some (ns,bh) 
systems should have been detected already, this estimate 
serves as a lower bound for effects that select against the 
detection of these binaries. 



2 ANALYSIS 

To estimate the frequency of (ns,bh) systems in the Galaxy, 
we construct a probability for a single binary system to be 
generated with the initial parameters that are a necessary 
(but not sufficient) condition for (ns,bh) formation. We as- 
sume that stars with masses between 8 and 25 Mq become 
neutron stars, with stars above this limit bec oming black 
holes. This is consisten t with recent work by iFrver et alJ 
j2002f) and lHeeer et al J ]2003f) on the end state of massive 
stars. The binarity fraction is taken to be 50% (i. e., 2 out 
of 3 stars are in a binary), and a standard Salpeter IMF is 
assumed (dN = m~ 2 ' S5 dm), drawn from 8 to 100 M . The 
mass of the secondary is determined by drawing the mass 
ratio q from a flat distribution. The choice of a mass ra- 
tio distribution is still quit e contentious; a study of massive 
stars in the Orion nebula (Preibi sch et alj H-999) showed a 
secondary IMF that is consistent with that of the field star 
population, but a flat distribution could not be ruled out 
either. Given the inherent obs ervational biases in studying 
systems with wide mass ratios llHogeveenlll992h . a flat mass 
ratio distribution is not an unreasonable assumption here. A 
distribution skewed towards lower mass secondaries would, 
of course, result in a smaller number of (ns,bh) systems. 

To form a (ns,bh) system via this channel, a number 
of conditions must occur in order, each with an associated 
probability. First, the system must have a small enough 
initial separation for a first mass transfer episode. Second, 
enough mass must be transferred from the primary to the 
secondary to ensure that the secondary becomes a black 
hole, with a reduction in its evolutionary timescale com- 
mensurate with the mass gained. The primary must retain 
enough mass to become a neutron star. We assume that 
the main sequence lifetime of the primary is not extended 
by mass loss (as it is the envelope which is lost). Third, 
the system must survive the natal kick that accompanies 
the collapse of the primary to a neutron star. Fourth, the 
perturbed system must retain a sufficiently close orbit that 
a second bout of mass transfer can occur, this time from 
the more massive secondary onto the neutron star primary. 
Lastly, the system must survive any natal kick arising from 
the collapse of the secondary into a black hole. 

Setting aside for a moment the question of whether 
many (ns,bh) progenitor systems survive the first natal kick, 
we can address in a basic, analytical way the potential upper 
limit for the formation rate of these binaries. 

The probability that an arbitrary binary will transfer 



mass as the primary leaves the main sequence depends upon 
the ZAMS mass Mo of the primary, the initial mass ratio go, 
and the initial semi- major axis ao- Practically speakin g, such 
syste ms have initial periods of less than 10 years iPols et alJ 
1991), with semi-major axes on the order of a few thousand 
solar radii. Initial orbital separations are typically tak en to 
be di stributed evenly in log a, out to roug hly lO 6 ^ HE3 
1983). As the region of interest occupies three decades in a, 
out of a possible six, we take the probability of initial mass 
transfer to be 1/2. 

We assume that the entire envelope of the primary is 
transferred during this initial episode, though based on the 
degree to which the transfer is non-conservative, not all of 
this matter will be accreted by the secondary. With a mass 
loss fraction /, the final masses (Mi, M2) of the components 
are related to their initial masses (Mi,M°) by the expres- 
sions 

Mi = M core (1) 
M 2 = M 2 ° + (M 1 ° - M core ){l - f) (2) 

where M core is the core mass of the primary, wh ich we take 
to be 0.058(M 1 °) 1 ' 57 M Q dlben fc Tutukovll 19851) . 

The mass loss fraction / is taken to be a linear 
function of the initi al mass ratio only JPols et alJ Il99ll : 
iPorteeies Zwartll995f) . ranging from 1 at go = 0.2 to when 
go = 0.6. Values above or below this range assume a value 
of unity or zero, respectively. An alternative is to treat / as 
a quadratic function of g; in the simplest case, / = 1 — g 2 . 
We will have more to say about the relevance of this param- 
eter later when discussing the results of a detailed binary 
evolution code in section |H] 

The final probability to consider is whether the system 
will remain sufficiently closely-bound to transfer mass from 
the black hole progenitor back onto the nascent neutron star. 
To first order, we consider that any system that survives 
the first natal kick will eventually go on to a second bout 
of mass transfer. Given the extreme difference between the 
imparted kick speed and typical orbital speeds for this sce- 
nario, those systems receiving a prograde kick will almost 
invariably be disrupted, while systems that remain intact 
may still be close enough for mass transfer after tidal circu- 
larisation has completed. 

The question of normalisation must also be consid- 
ered. We choose to normalise to the Galactic supernova rate 
of 0.01 yr -1 (consiste nt with the 0.008 yr _1 reported by 
ICappellaro et afll997l) . We assume all stars with a zero-age 
main sequence mass between 8 and 25 Mq contribute to the 
supernova rate. The question of whether black holes experi- 
ence a supernova at formation is a complex and unresolved 
one, with a critical dep endence on the star's pre-collapse 
mass iFrver et al] 12003) . For our purposes, we ignore the 
natal kick imparted to a nascent black hole, and do not 
count these events towards the supernova rate. 

From the above considerations, we can establish the 
range of primary and secondary masses for which the (ns,bh) 
scenario is physically possible. Figure shows the bound- 
aries in mass space for (ns,bh) candidates. The four bound- 
ary conditions are indicated by roman numerals, and repre- 
sent the following: (I) that the envelope mass M e of the pri- 
mary is greater than the difference between the ZAMS mass 
M2 of the secondary and 25 Mq (the threshold for black hole 
formation), so that M e = Mi - 0.058M1 1 ' 57 > 25 - M 2 . (II) 



Recycled pulsars with black hole companions 3 



The tautological requirement that the initial primary mass 
is greater than that of the secondary (i. e., < q < 1). 
(Ill) The ZAMS mass of the primary does not exceed the 
threshold for black hole formation (Mi < 25M Q ). (IV) If 
the initial mass transfer is very non-conservative, the enve- 
lope mass M e will not be sufficient to make the secondary 
a black hole, so we require that AI2 > 25 — M e (l — /). Note 
that this last requirement is essentially a more stringent ver- 
sion of condition (I) , and only comes into play for primaries 
above 20 M©. 

A crude estimate of the (ns,bh) formation rate can be 
obtained by simply calculating the area bounded by the 
four aforementioned conditions, weighted by the initial mass 
function. This gives an upper limit to the formation rate, as 
it ignores the effects of natal kicks, and the likelihood of 
the two mass transfer episodes. From our previous assump- 
tions (Salpeter IMF, flat mass ratio distribution, and a 50% 
binarity), and normalising to the Galactic SNR, we find a 
formation rate of 1.6 x 10 _4 yr _1 . If we consider a more 
severely nonconservative mass transfer, with / = I — q 2 , the 
formation rate shrinks by a factor of two, to 7.5 x 10 yr -1 . 



3 DISCUSSION 

Having determined the formation rate of recycled (ns,bh) 
binaries, we can estimate the number of such systems per 
pulsar by comparing the product of the respective formation 
rates and pulsar lifetimes. This is dependent upon the ex- 
tent to which the pulsar is recycled, a complicated topic 
which we do not address at present. We take a lifetime 
ti = 10 s yr as a reasonable average, consistent with re- 
cycled pulsar lifetimes inf erred from neutron star binaries 
iArzoumanian et alJll999l) . For unrecycled pulsars, we as- 
sume an average lifetime of 2 x 10 7 yr. The results are sum- 
marised in Tabic 

The estimate of the (ns,bh) formation rate above is 
clearly too large; among approximately 1500 known pulsars, 
many of these systems should have already been seen. Most 
of the putative (ns,bh) systems should in fact have been 
disrupted by symmetric and asymmetric kicks during the 
formation of the neutron star. To better understand the ex- 
tent to which kicks prevent these systems from forming, we 
employ a more detailed numerical model of binary evolution, 
including several choices regarding the distribution of natal 
kick speeds. 

We chose to use t he SeBa package for our numerica l 
model, as described in iPorteeies Zwart fc Verbunti J1996). 
However, several design choices require comment at this 
point. Our prescript ion for handling helium star winds is 
that of lLangerl £l989) . The stellar evolution tracks that form 
th e core of the pa c kage a re the analytic interpolations given 
in lEggleton et al.1 dl989l) . 

Probably the least well-constrained quantity in our 
model is the initial black hole mass function; i. e., the map- 
ping of pre-collapse masses and the mass of the resulting 
black hole. The approac h used in SeBa is based on that of 
iFrver fc Kalogeral 112001 ). There, the total binding energy 
of the progenitor at collapse is compared to the energy re- 
leased in the supernova, multiplied by an efficiency factor 
which is taken to be roughly 1/2. In SeBa, we take the 
supernova energy as fixed at 10 51 ergs, instead of a func- 



tion of the progenitor mass. The mass of the CO core is 
taken as the minimum mass of the black hole, to which 
some fallback material may be added. The binding ener- 
gies for the helium and hydrogen envelopes are calculated 
separately, and then each is compared to the supernova 
energy in turn. If the hydrogen envelope has a mass Mh 
with binding energy Eh , a supernova providing energy Esn 
will produce an amount of fallback hydrogen Fh given by 
F H = M H x (1 - (0 < F H < 1). Any remaining su- 

pernova energy after the hydrogen layer is unbound is then 
applied to the helium layer, so that the amount of fallback 
helium is just F He = M He x (1 - Es %~f H )- K any of the 
hydrogen falls back on to the black hole, the entire helium 
envelope will also be captured, by definition. The final black 
hole mass is then just the CO core mass plus the total fall- 
back material. 

To better investigate the role of kicks, we chose three 
separate kick distributions. The first allows a symmetric kick 
only, with a zero-magnitude asymmetric kick. Next, we per- 
formed ajenesof runs usin g the asymmetri c kick distribu- 
tion of lPaczvnskil Jl990T) and lHartmanl Jj.997), with a disper- 
sio n speed of 300 km s - 1 . Las t we use the model formulated 
bv lCordes fc Chernofil il998T) . derived from observations of 
the Galactic pulsar scale height. This last distribution is a 
double-gaussian, drawing 86% of the velocities from a gaus- 
sian with a dispersion of Vdisp = 175 km s _1 , and the re- 
mainder from a gaussian with Vd iap — 700 km s _1 . 

The results of this simulation are summarised in Table 
directly below the event rates derived from a direct inte- 
gration as described in the Analysis section above (sets A 
and B in the table) . The set C come from the output of the 
SeBa binary evolution code with no asymmetric kicks ap- 
plied at neutron star formation (symmetric kicks still occur). 
D and E also come from the SeBa routine, but with natal 
kicks drawn, respectively, from the Paczynski and Cordes- 
Chernoff distributions described above. The (ns,bh) forma- 
tion rate is found by scaling the output to the Galactic su- 
pernova rate of 0.01 yr -1 . As well, it is assumed that recycled 
pulsars will live an average of 10 8 yr, as before, while slow 
pulsars are taken to have a lifetime of 2 x 10 7 yr for this 
estimation. We immediately see the enormous impact that 
natal kicks have on the formation rates of such binaries. 
The considerable difference between the analytical results 
(A and B), and the zero-kick numerical simulations arises 
from ignoring the effects of wind mass loss in the former. 
Systems that survive the first supernova may widen con- 
siderably due to wind mass loss from the secondary as it 
leaves the main sequence, preventing a second mass trans- 
fer episode iPortegies Zwart et al.lll997l) . Hence our initial 
overestimate of the formation rate. 

Given that the number of known pulsars is less than 
two thousand, it is not at all surprising that no (ns,bh) sys- 
tem has been found to date. While strongly dependant upon 
the chosen natal kick distribution, it is clear that the total 
number of known pulsars must grow by a factor of several 
before even one such binary is likely to be found. 

Figure|2]shows the initial ZAMS masses of both compo- 
nents of (ns,bh) progenitors. For each of the three choices of 
kick distribution, 1.2 x 10 6 binaries were generated, with the 
remaining intact (ns,bh) systems shown on the plot. We then 
overlayed the criteria from Figure to check how well our 
constraints were borne out. The bulk of the (ns,bh) systems 
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fall inside the theoretical bounds we established. One excep- 
tion is our assumption about the extent of nonconservative 
mass transfer as a function of the mass ratio q. Clearly, our 
initial linear approximation for the conservation factor / was 
too lenient. The dashed line in Figures and [5] represents 
the mass conservation condition when / = 1 — q 2 . This as- 
sumption better matches the data, but the true curve lies 
between these extremes. Of special note are the scattering 
of (ns,bh) systems with very high primary ZAMS masses, 
well to the right of condition III in both figures. These rare 
systems do result in (ns,bh) binaries, but are the result of 
Case A mass transfer from a very massive primary to a sec- 
ondary that is just below the mass threshold for black hole 
formation. We did not consider Case A mass transfer when 
formulating our initial scenario, but it is interesting to note 
this alternate, albeit secondary formation channel. The frac- 
tion of (ns,bh) systems arising from Case A mass transfer 
varies as a function of the natal kick distribution. With no 
natal kicks, this channel comprises about 10% of the total 
(ns,bh) population. This rises to roughly 30% when Cordes- 
Chernoff or Paczynski-type kicks are applied, since systems 
that exhibit Case A mass transfer are closely-bound by def- 
inition, and more likely to survive a natal kick. 

It is also interesting to consider the "normal" formation 
channel of (bh,ns) binaries; i. e., where the black hole forms 
first. The final two columns of Table show the formation 
rate of these systems, and an estimate of the ratio of regular 
pulsars to (bh,ns) binaries. While such systems form roughly 
two orders of magnitude more frequently than (ns,bh) bina- 
ries, the shorter lifetime of the unrecycled pulsar means that 
(bh,ns) are only about 20-30 times more common than the 
reversed channel. Looking at the expectation value of pul- 
sars to (bh,ns) systems, we note that the number for set E is 
comparable to the current number of known pulsars, imply- 
ing that detection of such a system could occur in the very 
near future. 
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Figure 1. Shows the range of stellar mass pairs that could con- 
ceivably lead to a (ns,bh) end system. The four boundary condi- 
tions are delineated by Roman numerals, and are derived from: 
(I) the requirement that a sufficient amount of mass be trans- 
ferred from primary to secondary to ensure that the latter be- 
comes a black hole, (II) the tautological requirement that the 
primary ZAMS mass is greater than that of the secondary, (III) 
the maximum mass of the primary, above which the star must 
inevitably become a black hole (25 Mq), and (IV) the limit of 
non-conservative mass transfer at low values of q. As mass trans- 
fer becomes less conservative, a larger envelope mass is required to 
raise the secondary above the black hole mass threshold. Condi- 
tion (IV) results from an alternative assumption of more severely 
nonconservative mass transfer. 



4 CONCLUSIONS 

We investigated a channel for the formation of (ns,bh) bina- 
ries, where a bout of mass transfer reverses the ordinary evo- 
lutionary order, resulting in a neutron star which forms first. 
This neutron star can then be spun up by a second episode 
of mass transfer from the black hole progenitor star, result- 
ing in a recycled pulsar orbiting a black hole. The rapidly- 
spinning pulsar should live several times longer than a typ- 
ical pulsar. We then estimated the formation rate of such 
systems via both a simple analytical calculation and then a 
more detailed binary evolution code, and used assumptions 
about the relative lifetimes of the recycled pulsars to calcu- 
late the expected ratio of field pulsars to (ns,bh) systems. 
The result is that the current known pulsar population must 
grow by factors of several before such a system is likely to 
be found, despite containing a longer-lived recycled pulsar. 

Using the same methods, we show that normal (bh,ns) 
binaries should be more common than (ns,bh) systems by an 
order of magnitude, and that the expected number of pul- 
sars per (bh,ns) system is comparable to the current known 
pulsar population. For both (ns,bh) and (bh,ns) binaries, the 
strongest single factor in determining the formation rate is 



the distribution from which natal asymmetric kicks are se- 
lected. 
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Table 1. Summary of recycled (ns,bh) formation rates 



Set 


Kick Disp.AV 


(ns,bh) yr 1 


(ns,bh) Total 


PSR/(ns,bh) 


(bh,ns) yr 1 


PSR/(bh,ns) 


A 





1.6 x 10~ 4 


16000 


13 


6.5 X 10" 4 


15 


B 





7.5 x 10~ 5 


7500 


27 


6.5 X 10~ 4 


15 


C 





9.3 x 10" 7 


93 


2200 


1.4 x lO" 4 


70 


D 


300 


1.2 x 10~ 7 


12 


17000 


2.1 x 10~ 5 


470 


E 


175/700 


5 x 10~ 8 


5 


40000 


6.5 x 10~ 6 


1500 



Set A is the result of direct integration of the area in figure weighted by the IMF, with the mass loss via 
nonconservative mass transfer given by a simple linear relation, where B is the same calculation, with the mass 
loss factor varying as the square of the mass ratio q; sets C, D and E arise from detailed binary simulations using 
the SeBa package. Set C includes symmetric Blaa uw kicks but n o asymmetric kicks. Set D allows for a random 
asymmetric kick drawn from th e distribution of iPaczynskj ll99Cfl . with a dispersion of 300 km s^ 1 ; set E draws 
its kick from the distribution of lCordes Kr, ChernofTi^99^^ a double-Gaussian with dispersions of 175 and 700 km 
s — , respectively. The (ns,bh) formation rate gives the number of recycled (ns,bh) systems per year, while the next 
column estimates the total number of recycled (ns,bh) systems in our Galaxy, by multiplying the formation rate by 
the lifetime of a typical recycled pulsar, taken here to be 10 s yr. Column 5 shows the expected number of pulsars 
per single (ns,bh) system. Here we assume that pulsars are formed at the supernova rate, and live an average of 
2 x 10 7 years without recycling. The final two columns show the formation rate of (bh,ns) systems, and the number 
of pulsars expected for every (bh,ns) system. 
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Figure 2. Primary and secondary ZAMS masses for recycled 
(ns,bh) systems generated by the SeBa stellar evolution code. 
Dots represent systems formed when no asymmetric kicks apply; 
triangles and squares show systems that survived a Paczynski- or 
Cordes-Chernoff-type natal kick, respectively. These systems are 
principally confined to the region bounded by the four theoretical 
conditions. Outliers at low masses result from a lower mass limit 
for black hole formation. Several confirmed instance of Case A 
mass transfer also resulted in the formation of recycled (ns,bh) 
binaries (the high primary mass outliers), though our putative 
scenario initially considered only Case B mass transfer for this 
channel. 
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